Measurement of tt Spin Correlation in pp Collisions Using the CDF II 

Detector at the Tevatron 
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The ti spin correlation at production is a fundamental prediction of QCD and a potentially incisive 
test of new physics coupled to top quarks. We measure the tt spin state in pp collisions at y^i = 1.96 
TeV using 1001 candidate events in the lepton plus jets decay channel reconstructed in the CDF II 
detector. In the helicity basis, for a top-quark mass of 172.5 GeV/c^, we find a spin correlation 
coefficient k = 0.60 ± 0.50 (stat) ± 0.16 (syst), consistent with the QCD prediction, k ~ 0.40. 

PACS numbers; 12.38.Qk,14.65.Ha 



In quark-pair production by the strong inter- 
action, the quark spins are entangled according 
to the short distance dynamics of quantum chro- 
modynamics (QCD) The spin state is observ- 
able in angular correlations among the quark de- 
cay products induced by the V-A (Vector minus 
Axial- vector) nature of the weak interaction but 
is typically destroyed by the depolarizing effects 
of hadronization before the decay can proceed. 
The top quark is an exception to this rule. Be- 
cause of its large mass, the top-quark lifetime is 
shorter than the fragmentation timescale, cutting 
off the long distance QCD effects and transmit- 
ting the tt production configuration to the final 
state. Measurement of the tt spin configuration 
is a first look at a bare-quark pair at produc- 
tion. The measurement tests the fundamental 
predictions of QCD [ll-IH and could be a sensi- 
tive discriminant of new physics coupled to top 
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quarks 0, For example, a tt resonance ap- 
pearing as an excess in the invariant-mass spec- 
trum can be verified as a Kaluza-Klein graviton 
through measurement of the spin correlation as 
described in Ref. [7|. 

Because final state charged leptons have the 
strongest correlation to the top-quark spin, the 
tt spin correlation is usually discussed in terms 
of the dilepton final state tt ^ {W+b){W-b) 
{£i'){£'i'')bb [4]. This mode suffers from a small 
branching ratio and poor definition of the top- 
quark kinematics due to the presence of two un- 
detectable neutrinos. A previous measurement 
of the tt spin correlation was limited to a small 
sample of just six events in this mode 

We report on a new measurement of the tt spin 
correlation in pp collisions at 1.96 TeV with a 
data sample corresponding to an integrated lu- 
minosity of 4.3 fb"'^ collected with the CDF II 
detector at the Fermilab Tevatron. We mea- 
sure the spin correlation of pair-produced quarks 
for the first time in the lepton plus jets decay 
topology, tt -J> {W+b){W-b) _^ {udb){ivb) or 
tt {W+b){W-b) {£vb){udb) [§]. In this de- 
cay mode, we take advantage of a large branching 
ratio compared to the dilepton channel and the 
well-constrained ti kinematics in the lepton plus 
jets final state with only one neutrino. The mea- 
surement relies critically on a new technique for 
identifying the final state down-type quark {d or 
s), which has the same spin-analyzing power as 
a charged lepton. We expect the spin correla- 
tion measurement to show the dominance of ti 
production via the J=l qq annihilation channel 
that occurs in ~85% of pp collisions at the Teva- 
tron [13. 

We work in the helicity basis, where the spin- 
quantization axis is defined as the direction of 
motion of the t (or i) quark in the ti rest frame. 
There are other quantization axes which predict 
a larger value for the spin correlation [3|, but 
they do not provide any significant increase in 
the statistical sensitivity of our approach, so we 
work with the simpler helicity basis. A quark 
is called right-handed (ti^)/left-handed (t^) if its 
spin is oriented along/opposite to its direction 
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of motion. In the tt rest frame the quarks move 
back-to-back; thus the same-spin states with J=l 
are those with opposite hehcity: ix^fl and tf/tL. 
Near the energy threshold for tt production, the 
opposite-helicity fraction is predicted in the stan- 
dard model (SM) to be ~67% for tt production 
via qq annihilation, while for top quarks with 
large momenta compared to the top-quark mass, 
helicity is approximately conserved and this frac- 
tion rises to ~100% Integrating over all 
top-quark momenta according to the parton dis- 
tribution functions and adding the small ('^15%) 
J = contribution from gluon-gluon fusion pro- 
cesses, we expect to find an opposite-helicity frac- 
tion [l|, n 



FoH — 



<7{tRtL) + Cr{tLtR) 



<j{tRtR) + a{tLtL) + (j{tRtL) + (j{tLtR) 

0.70. (1) 



Fqh is simply related to the spin correlation 
coefficient k that measures the fractional differ- 
ence between the number of events in which the 
top-quark spins are aligned and the number of 
events in which they have opposite directions: 
K = 2FoH - 1. We thus expect k « 0.40 P, H, 
while for uncorrelated spins, k — 0.0 and Fqh = 
0.5. 

In top-quark decays in the SM the V-A cou- 
plings fix the angular distributions of the decay 
products according to the polarization of the par- 
ent top quark via 



1 dV 



r d cos i 



-(1 ± A^cos^i), 



(2) 



where the positive/negative sign is used for right- 
handed/left-handed quarks, and the helicity an- 
gle 9i is defined as the angle between the spin- 
quantization direction and the momentum of the 
decay particle in the rest frame of its parent 
top quark. In the V-A weak decay, the spin- 
analyzing-power coefficient Ai is equal to -1-1.0 
for the charged lepton or down-type quark, -0.41 
for the bottom quark, and -0.31 for the neu- 
trino or up-type quark, with the signs reversed 
for antitop-quark decays Q . The ti spin correla- 
tion connects the daughter helicity angles on each 
side of the decay. The differential cross-section in 
these variables is 



KAiAj cos 9i cos6j 



a d{cos9i)d{cos 6j) 



(3) 

where i and j refer to top-quark and antitop- 
quark decay products respectively [3l|. 

For each of the four possible tt helicity states, 
we create model templates for the distributions 
of cos 9 1 cos 9d and cos 9i cos 9^, where the charged 
lepton I is a decay product from one top quark 
in the pair and the quarks d and b are decay 
products from the other quark. We then find the 
relative normalization of these model templates 
that gives the best fit to a two-dimensional distri- 
bution of these variables in the data. The model 
templates account for all acceptance effects and 
dilutions due to event reconstruction, so that the 
parton-level value of Fqh follows directly from 
the template fit to the data. 

CDF II [ll| is a general purpose, az- 
imuthally and forward-backward symmetric de- 
tector. Charged-particle directions and momenta 
are measured with a silicon tracker [Hj and a 
drift chamber [l^ in a 1.4 T solenoidal magnetic 
field. Electromagnetic and hadronic calorime- 
ters [31 are located beyond the solenoid and al- 
low for jet and missing Et reconstruction. Be- 
yond the calorimeter, muon chambers [l5| pro- 
vide coverage for the pseudorapidity range \r]\ < 
1.0. We use a cylindrical coordinate system with 
its origin at the center of the detector and the z 
axis along the proton direction [l6j . 

Lepton plus jets events are selected by requir- 
ing one electron or muon with transverse momen- 
tum of at least 20 GeV/c and \ri\ < 1.0, missing 
transverse energy of at least 20 GeV, and four or 
more jets with transverse energy of at least 20 
GeV and < 2.0, at least one of which must be 
tagged as a 6 jet by the presence of a displaced 
secondary vertex [l3|. This selection yields 1001 
total candidate events, 224 of which have two 
tagged b jets. 

Non-tt backgrounds are well-constrained by 
precision tt cross-section measurements 18| , with 
a predicted total of 215 ± 48 background events. 
Non-tf models are checked against background- 
enriched sidebands with no tagged b jets and are 
found to give very good representations of the 
normalizations and kinematics in all variables, 
including lepton and jet energies and angular dis- 
tributions. 

The helicity angles are determined in a com- 
plete reconstruction of the tt kinematics in tt ^■ 
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{Wb){Wb) — 7> {ivh){udb), where we constrain 
M{tv) = M{ud) = 80.4 GeV/c^, the mass of 
the W boson, and M{ivb) = M{udb) = 172.5 
GeV/c^, the top-quark mass, and require any 
tagged b jets to be identified with b partons. The 
constraints were chosen to be close to the world 
averages in Ref. fl9|. Each of the 24 possible 
jet-to-parton assignments is evaluated using a 
comparison to the tt hypothesis with the above 
constraints, and we choose the assignment with 
the lowest value [l^l . This procedure correctly 
assigns all jets to the corresponding partons in 
approximately 37% of events. All effects of an- 
gular acceptance and jet reconstruction and mis- 
assignment are fully modeled by our simulated 
samples. 

Down-type-quark identification relies on the 
V-A decay correlation that tends to send the 
down-type quark in the direction opposite that 
of the hadronically decaying W boson in the top- 
quark rest frame. We therefore assign the down- 
type quark as the jet that, in the W boson rest 
frame, is closest to the bottom jet identified as 
coming from the same top quark as the W bo- 
son [3|. Simulation studies show that this al- 
gorithm correctly identifies the down-type quark 
60% of the time. 

The same-helicity and opposite-helicity model 
templates are created with a customized version 
of the HERWIG event generation software pack- 
age dlj that implements the angular distribu- 
tion of Eq. [5] for the charged lepton or down- 
type quark, with a tunable choice of right- or 
left-handed top quarks, and preserves all the 
other expected spin correlations [2^. We cre- 
ate four different simulated samples, correspond- 
ing to the four possible top-quark-pair helic- 
ity states: tLtn, intL, hth, and tRtR. QCD 
interactions respect both the parity symmetry 
(P) and the combined symmetry of parity and 
charge conjugation (CP). Because CP trans- 
forms tflifl — ?> iL^L, we can define the same- 
helicity (SH) model template shape to be the 
symmetric sum of a{tjitji) -f (T(ii,ii,). Since P 
transforms iutL ihtR, we let the opposite- 
helicity (OH) model template shape be the sym- 
metric sum of cr(tiitL) + uitLtft). 

Figure [T] compares the SH and OH model tem- 
plates after detector simulation, event selection, 
and reconstruction in the two distributions that 
we use for the measurement, cos 9i cos 0d and 
cos^i cosdb- Our sensitivity results from the SH 
model template being shifted towards negative 
values of cos 0; cos 0^, while the OH model tem- 



plate is shifted towards positive values, with the 
opposite shifts occurring in the cos 6i cos 9b dis- 
tribution. 
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FIG. 1: Distributions of the cos 6; cos and 
cos^; cos Ob variables, after detector simulation, event 
selection, and reconstruction, in our same-helicity 
and opposite-helicity simulated ti samples. 

We perform our measurement using a binned 
likelihood fit to find the relative normalization 
of these model templates that gives the best 
simultaneous representation of cos 9i cos 9d and 
cos 9i cos 9h in our data. The background normal- 
ization is constrained to be close to the predicted 
value, with a Gaussian uncertainty, but the same- 
helicity fraction Fsh and opposite-helicity frac- 
tion FoH are allowed to float freely. We do not 
require that Fsh and Fqh be constrained to 
physical values between and 1, but we do re- 
quire FsH+FoH = 1- The fit runs over all bins in 
a two-dimensional distribution of cos 9i cos 9d vs 
cos 9 1 cos 9b. The expected statistical uncertainty 
for Fqh is approximately 0.23, corresponding to 
an uncertainty for n of 0.46, and is independent 
of the actual value of Fqh and k. 

Additional contributions to the uncertainty 
result from incomplete knowledge of the back- 



6 



TABLE I: Systematic Uncertainties on Foh 



bystcmatic 


Uncertainty 


Generator dependence 


0.060 


JES 


0.042 


ISR/FSR 


0.030 


Background shape 


0.023 


Color rcconncction 


0.009 


PDF 


0.007 


Parton shower 


0.006 


Background size 


0.002 


Total uncertainty 


0.083 



ground size and shape, of the exact detector 
response, and of the parton distribution func- 
tions (PDF), and are estimated by performing 
the measurement in simulated samples with rea- 
sonable variations in the model assumptions. 
These systematic uncertainties are shown in Ta- 
ble m The largest uncertainty, generator de- 
pendence, results from small biases seen when 
testing with simulated samples created using a 
range of event generation software packages, in- 
cluding HERWiG '21], pythia [23|, alpgen [23 |. 
and MADEVENT [25]. Other significant contribu- 
tions come from the uncertainty of the jet energy 
scale (JES) during event reconstruction and un- 
certainty in the amount of initial and final state 
radiation (ISR/FSR) in our observed tt events. 
The small variation of Fqh with the assumed 
value of the top-quark mass is not included in 
our systematic uncertainty; our measurement as- 
sumes a mass of 172.5 GeV/c^ for the top quark. 

The final result of our fit to the two- 
dimensional distribution cos 6i cos 6d vs. 
cos 0; cos 6*6 is shown in Fig. [21 This figure 
shows one-dimensional distributions of both 
variables, with our data being compared to the 
sum of the background model, same-helicity 
model, and opposite-helicity model, with the 
model normalizations determined by our fit 
result. Assuming the top-quark mass is 172.5 
GeV/c^, we find an opposite-helicity fraction of 



Foh = 0.80 ± 0.25 (stat) ± 0.08 (syst) . 

Converting this to the spin correlation coefiicient, 
using K = 2FoH — 1, yields 



I I OH Model 
^ ■ SH Model 
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• Data 




-0.4 -0.2 0.2 0.4 
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K=0.60±0.50±0.16 
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Distribution of the cos di cos 9d and 
't, variables in data compared to the sum of 



FIG. 2: 

cos 9i cos 

our background model, the same-helicity model tem- 
plate, and the opposite-helicity model template. The 
relative normalizations of the model distributions are 
determined by our fit result. 



K = 0.60 ±0.50 (stat) ±0.16 (syst) . 

This first measurement of the top quark-pair 
spin correlation in the lepton plus jets decay 
channel agrees well with the theoretical predic- 
tion of K w 0.40 [J, 3, although the statistical 
uncertainty is still large. Simulated experiments 
with larger datasets indicate that if the Teva- 
tron dataset reaches 15 fb^^ before the end of 
the Tevatron lifetime, the expected statistical un- 
certainty on K would be reduced to 0.26. This 
technique can thus be applied in future measure- 
ments with larger datasets collected at the Teva- 
tron and LHC to constrain the production spin 
structure or to connect with other anomalies that 
may show up in the reconstructable ti kinematics 
of the lepton plus jet sample. 
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